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Abstract 

We consider Wong equations for a particle with a continuous mass 
spectrum in a random Yang-Mills field approximating the quantum field 
at finite temperature. We show that particle time evolution can be ap- 
proximated by a relativistic diffusion. Kubo's generator of the relativistic 
diffusion is defined as an expectation value of the square of the Liouville 
operator. We discuss hydrodynamics, energy and entropy of the gas of 
diffusing particles. 

1 Introduction 

In non-relativistic physics the dynamics of a test (tracer) particle in a medium of 
some other particles is often approximated by diffusion [1]. We encounter such 
a problem in high-energy physics as well. The diffusion approximation has been 
applied to a description of the quark-gluon plasma and to heavy ion collisions 
[2] [3] [4] [5] . If QCD is considered as the fundamental theory of strong interactions 
then the diffusion should result from the QCD transport theory [6] [7] [8]. The 
form of the diffusion operator can in principle be determined on the basis of the 
Boltzmann kinetic equation after a calculation of the scattering amplitude [3] [9] . 
A semi-classical approach to the transport theory based on the Wong equation ( 
derived from quantum dynamics in [10] [11] [12]) leads to an analog of the Vlasov 
equation [6] [7] [11] [13] . Quantized gauge fields in the Vlasov equation have been 
discussed in [14] [15] [16]. These authors have shown that quark's time evolution 
in a quantum gauge field leads to a diffusion of momentum and color. However, 
the generator of the diffusion has not been determined in an unambiguous way. 
It is well-known [17] [18] that the Vlasov equation describing a tracer particle 
moving in plasma in a random electromagnetic field produced by a chaotic 
motion of other particles can be approximated by a diffusion equation. The 
generator of the diffusion has been calculated in [17] [18]. A detailed discussion 
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of the diffusion approximation to a random non-relativistic dynamics can be 
found in [19] [20] and [21]. 

The mathematical theory of a relativistic diffusion is not so well developed 
as the non-relativistic one (see the reviews in [22] [23]). A relativistic diffusion in 
the configuration space does not exist [24] [25] . If we assume that the diffusing 
particle has a fixed mass then an analog of the Kramers diffusion is uniquely 
defined by its generator being the Laplacc-Bcltrami operator on the mass-shell 
as shown by Schay [26] and Dudley [27] . There are some other suggestions for a 
definition of the relativistic diffusion without the assumption that the diffusion 
takes place on the mass-shell [22] [23] [28]. 

In this paper we start from the basic principles of quantum gauge theory at 
finite temperature. We assume the gauge and Lorentz covariance in the most 
general form. It is known that the thermodynamic equilibrium is defined in 
the rest frame. Then, if the finite temperature theory is Lorentz covariant then 
it can be transformed to an arbitrary Lorentz frame described by a unit time- 
like four- vector . We assume that the two-point function of gauge invariant 
field strength [29] [30] F^ u (x) depends on two four- vectors the position x and 
the frame w. It comes out that the diffusion generator depends only on the 
form of the two-point function at coinciding points. The diffusion tensor can 
be related to the energy-momentum tensor of gluons in the thermal state. We 
do not assume that the particles (quarks) are on the mass-shell as we did in 
the case of (zero-temperature) QED in [31]. Only free stable particles satisfy 
the mass-shell condition. Relaxing this condition allows an explicit use of the 
Lorentz covariance ( four independent components of the momcntum;the trans- 
port theory for particles with a continuous mass spectrum has been discussed 
also in ref.[32]). The plan of the paper is the following. In sec. 2 we formulate 
the Wong equation in the random field. In sec. 3 we calculate the expectation 
value of the random evolution (the square of the Liouville operator) which de- 
termines the diffusion generator. In sec. 4 we take the diffusion model resulting 
from the quantum gauge theory at finite temperature as the starting point of 
the statistical mechanics of diffusing particles. We define the energy, free energy 
and entropy satisfying the second law of thermodynamics. We initiate a hydro- 
dynamic description of the gas of diffusing particles. The aim is to develop a 
macroscopic hydrodynamics on the basis of quantum gauge theory to supple- 
ment the phenomenological approach of Israel and Stuart [33] usually applied 
in heavy ion collisions [34]. 



2 Wong equation in a random Yang-Mills field 

The dynamics of a particle in a Yang-Mills field A® is described by the equations 
[10] 



dr ' 



f abc p»AlQ c , (3) 



where fj, = 0, 1, 2, 3 . It follows from eq.(2) that 

^«^)=0. (4) 

where »7 M „ = (1, — 1, — 1, —1). From eqs.(l) and (4) it follows that r is the proper 
time 

dr 2 =dx^dx^p- 2 . (5) 

If we consider an observable as a function (p on the product of the phase 
space and the internal (color) space (x,p) x SU(n) then while the coordinates, 
momenta and charges (x(t),p(t),Q(t)) evolve in the proper time according to 
eqs.(l)-(3) a function of observables <f) T = <P(x(t) , p(r) , Q{t)) evolves as 

d T 4> =P*$r + QaF^p^ - rV^^f- (6) 

In eq.(6) we assumed that po is an independent variable. This means a con- 
tinuous mass spectrum as in a description of unstable particles (resonances) in 
quantum field theory [32]. The Wigner function approach to the kinetic equa- 
tion [9] [35] also leads to the continuous mass spectrum well-known from the 
Kallen-Lehman representation of expectation values of quantum fields. We de- 
fine an expectation value of the observable in a state (probability distribution) 
ft as 

ft(</>) = (ft, <j>) = [ dxda(p)Q4>. (7) 



where the momentum distribution a will be expressed in the form (we assume 
m 2 > 0) 

da{p) = dpdm 2 p( m 2 )5(p 2 - m 2 c 2 ) (8) 

In order to describe a particle with a fixed mass m it is sufficient to choose 
p(m 2 ) = S(m 2 — Too). 

The evolution of ft is defined as an adjoint of the evolution of cj> 

(ft r ,0) = (ft,0 T ). (9) 

Then 

d T n T = -p^ - Q a F« uP »^ + f^Q aP »A b ^. (10) 

The requirement of the r- independence of the probability distribution (9 T ft T = 
0) gives the kinetic equation in the laboratory time t. This is the same equation 
as the one which can be derived by an elimination of r in favor of t in the 
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evolution equations (l)-(3). It will be useful to work with a matrix notation. 
Let T be a representation of the algebra of the group SU(n) 

[T a ,T b ] = f abc T c 

normalized by Tr(T a T b ) = -5 ab . We define Q = Q a T a and A p = T a A^. We 
consider the path-ordered exponential along a curve 7 

E/ 73 , = r(cxp(- J ' A^drf)). (11) 

with y as the endpoint of the curve. Then, 

-^U iy = -A,(y)U Jv . (12) 

The field strength F is constructed from A^. This assumption is expressed as 
the Bianchi identities 

e^D„F ap = 0, 

where D p F ap = d p F ap + \An,F„ p ]. We consider random gluon fields as an 
approximation to quantum fields. In general, correlation functions of gauge 
fields being gauge dependent are not well-defined. We are going to define a 
two-point correlation function of gauge invariant variables 

F lxu { 1 ,y) = U-^F llv {y)U iy . (13) 

F satisfies the Bianchi identities with an ordinary derivative (instead of the 
covariant one;this has been shown first in [29], see also [30]) 

e^Pd p F ap = (14) 

The Yang-Mills equations are [30] 

8' M F IJlv { 1 ,y) = U-^J v {y)U lv . (15) 

In terms of the path ordered exponential the solution of eq.(3) is 

Q(t) = U p (t)Q(0)U- 1 (t), (16) 

where U p (t) — U ly is calculated along the straight line in the direction p,i.e., 
7(s) = x — ps with y = x — pr . 

We define the corresponding gauge invariant field strength transported along 
the straight line 

F p „(p,x-ps) = U p {s)- 1 F pu {x-ps)U p {s). (17) 

We treat the Wong equations (l)-(3) in a random Yang-Mills field as an 
approximate description of the dynamics of quarks in quantum gauge theory at 
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finite temperature. Then, the expectation value of an observable O is defined 

(0)p = Tr(exp{-pw ll P' 1 )o) (18) 

where is the generator of space-time translations and the four-vector 
(w^w^ — 1) describes the reference frame [36]. In QED (and perturbative 
QCD) the expectation value (18) leads to the distribution of photons (gluons) 
according to the covariant Bose-Einstein law 

n(p;w) = (exp(0w v p v )-l)- 1 (19) 

The two-point correlation functions of gauge invariant variables F(p, x — ps) are 
invariant under transformations of the Poincare group (acting upon x,p and to) 
and under transformations of the (internal) SU(n) symmetry. We assume that 
the non-commutativity of quantum gauge fields at finite temperature can be 
neglected in the calculation of the particle time evolution. Then, the symmetry 
properties with respect to an exchange of indices together with the Bianchi 
identities (14) lead to the representation (see a more detailed discussion of the 
Abelian case in [31]) 



(F- V (p,x -ps))F b ap (p,x> -ps')) = S ab G^. ap (x - x'-p(s - s');s,s',w), 

(20) 

where 

G%- ap {y;s,s' ,w) 

= j dkG(s,s\k,p,w)(r]p a kvkp-rip P kvkp + riv P k cr kp-riv a kpkp)ex.p(iky) 

(21) 

G is a Lorentz invariant function of k, p and w. Let us note that 

G jL ;CT p(0; 0, °> w ) = VnvTvp - VnpTvp + ^ivpTap - V^Tpp (22) 

where 

Tpv = J dkG(0, 0, fc, 0, w)kp,k v = pw p w v - n E (rj^ - w M w„) (23) 

Eq.(23) gives a decomposition of the tensor T p v in terms of f]^ and w^Wu. if 
G is non-negative then p and tte must be non-negative. G at s = s' = does 
not depend on p because p can enter F only through the path starting in 
and ending in sp. T^ u is related to the expectation value (22) of the Yang-Mills 
energy-momentum tensor. In order to see this let us calculate 

—Tr{F tip F va r) pa - \^F ap F^) = 2T^ - \^ V TP (24) 

The expectation value of F pv can be different from zero if there are some source 
terms in the Lagrangian. So, adding the term jpA^ to the Lagrangian of QED 
gives 

{Fp,{x)) f 3 = j dyf(y)(d;G pv (x - y) - d x v G P p{x - y)) (25) 



5 



where G M „ is the two-point function of the potential A^. A gauge and Poincare 
invariant choice j M = in eq.(25) leads to 



= (27r) -3 / dk5(k 2 )5(k)n(k, w)(k a k„w^ - k^k^w^) = X(r]^ a w v - ri vcr w^) 



with a certain constant A ( rhs of eq.(26) satisfies the Bianchi identities (14)). 
In non-abelian gauge theories we may insert in the formula for the expectation 
value the Wilson loop (11) with the curve 7 M (s) = sw^ (this is the Polyakov loop 
applied for a description of confinement in finite temperature gauge theories) 



In a perturbative QCD we again derive eq.(26) with a certain constant A. 

In finite temperature QED p = w^T^w" in eq.(23) has the meaning of the 
photon energy density when G ~ 5(k 2 )n{k 1 w) . tte has the meaning of the 
pressure of photons. Such an interpretation is expected to hold true for gluons. 

3 Time evolution in a random gluon field 

In this section we study functions of particle trajectories. We consider kinetic 
evolution equations (6) which are of the form 



(26) 




(27) 



d s (f, = (X + Z + Y)<f> =(K + Y)t(>. 



(28) 



where K = X + Z. Let 



Y(s) = cxp(-sK)YcMsK). 



(29) 



Then, the solution of eq.(28) can be expressed as 



4> t = expitK)^ , 



(30) 



and 



a # = r(*)#. 

In eq.(28) we have according to eq.(7) 



(31) 




(32) 



Y = -Tr(QF^(x))p»^ 
Z = -f ahc Q a p»Al-£-. 



(33) 



(34) 



G 



Let us calculate exp(rK) (with K = X + Z). We have 

exp(rX)^ = exp(rX)^, (35) 

where 

drill = Z(T)tf. 

with 

Z(r) = C M-rX)ZcMrX) = -f^A^x-^Qa^. (36) 
From eq.(28) and (30) we obtain 

Y(s) = -Tr(Q(s)F, v (x + 
= -Tr(Q(0)F„ v (p,x sp))p^^- + s^-). 

We are going to approximate the evolution (28) by a diffusion. In general, this 
will be a diffusion on the product of the phase space and the group SU{n). 
As follows from eqs.(l)-(3) such an evolution is covariant but not invariant 
with respect to the gauge transformations. We show that the diffusion in the 
momentum space can be defined in a gauge invariant way. 

We apply the approach of Kubo [19]- [20] which approximates the random 
Liouville operator on the rhs of eq.(6) by the expectation values of Y and its 
square 

(#) ~ cxp (£ds(Y(s)) + l -{{ j\ 3 Y{ S )Y)) (38) 

where 

Y = Y — (Y) (39) 

If the mean value of Y(s) is different from zero then in the expectation value 
(<t>t) we could obtain the term of the first order in t 

K x = (Tr(QF^))p»d v (40) 

(from now on we denote the momentum derivatives by d v and space-time deriva- 
tives by ) However, such a term must vanish owing to the Poincare covariance 
because it is impossible to express the antisymmetric tensor (Tr(QF tlu )) by rf p 
and w M . Then, we have the second order term K 2 t 2 coming from the expansion 
of F^ v (x - sp) in s 

K 2 = {Tr{Qd x a F^))p»p°d v = \{r,^w v - r ]vr7 w^p"d v = XP^vTd* (41) 
where 

p/u, = ^ _ p -2 p n p . 
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The second order term resulting from the square of Y reads 



lim^o 4, ((/; Y(s)ds) 2 ) = lim T ^ & f Q ds f Q ds'(Y(s)Y (s>) + Y(s')Y(s)) 



Elementary calculations using eq.(22) give 

A w = -2K E q 2 d^p 2 d v + q 2 (p + TT E )((wp) 2 d^ - 2wpw^pPd ll d p 
+p 2 w tl w p d tl d p - w 2 p l, d v — 2wpw u d l/ ) 



(42) 



(43) 



where q 2 = Qa(T)Qa(i~) = Q a (0)Qa(0) is the square of the charge (which is a 
constant). The term sd^ of eq.(37) does not contribute to „4„,.This is so because 
it is of higher (third) order in r. The limit in eq.(42) depends only on DG(0) 
corresponding to s = s' — because taking into account s in G would give 
higher orders in r. As a consequence the coefficients do not depend on 
p because the dependence on p is always multiplied by s. According to Kubo 
[19]- [20] the t 2 behaviour at short time (small in comparison to the correlation 
time) of the random evolution goes over to the linear behaviour at large time 
with the same generator A w + K 2 . In order to justify a replacement of the 
Liouville evolution (6) by the diffusion (42) at an arbitrarily large time we have 
to apply the Markov approximation. In such a case the dynamics is viewed as 
a composition of independent short time evolutions ( see the discussion in Kubo 
[19]- [20] ). The rigorous proof [21] requires a proper scaling of fields and the 
time in order to determine precisely the conditions required for the limiting 
diffusive behavior. 

In the lowest order in the expectation value of the expansion of the solution 
(28) (besides the term (42)) we shall also have 

f ds f ds'{Z(s)Z(s')) (44) 
Jo Jo 

and 

[•T r-S 

ds I ds'(Z(s)Y(s')). (45) 



The term (44) would describe a charge diffusion in a Markovian approximation 
for the (AA) correlations. We would have to determine a form of the (AA) 
correlation function of the gauge field which would not be gauge invariant. 
The whole (non-perturbative) evolution exp(r-ftT) is gauge dependent although 
different gauges just rotate the charges (as can be seen from eq.(16)). In this 
sense fixing the gauge and approximating the charge evolution by the Brownian 
motion on SU(n) (a well-studied problem in the diffusion theory [37]) could still 
make sense (as discussed in [14] ) . The term (45) describes a mixing of the 
momentum diffusion with the color diffusion. In this paper we do not go into 
details of the color diffusion because of its complicated gauge dependence. 
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4 Relativistic diffusion at finite temperature 

The calculations of sec. 3 lead to the following diffusion equation in the proper 
time describing an evolution of a particle observable <j> ( a function of particle's 
position and momentum) 

d T <j> = g w <t> = P Kd*<j> - ad^p^ P 2 d^ - \ P 2 p^w u d^ 

+b({ P w) 2 d^d^ - (wp){w^pP + wPp")d p d^+p 2 w^wPd p d^) ( 46 ) 
-w 2 p v 8 v <j) - 2wpw v 8^(/) = p»d%<f) + A w <p 

where a = 2q 2 iTE , b = q 2 (p + TTE) and A is determined by the expectation value 
(41). 

The probability distribution evolves according to the adjoint equation (9). 
The independence of the proper-time parametrization {d T VL T = 0) is equivalent 
to the replacement of the proper time by the physical time t in the kinetic 
equation (10). This requirement gives the transport equation 

G+n = (47) 

where the adjoint is either in L 2 (dxda) or L 2 (dxdp)(it does not depend on a). 
Explicitly, 

pPdffi = d p (D p + \p 2 P^w u )tt = 8^- aP pv p 2 d v + \p 2 P pv w v + b{{pw) 2 d p 
-(wp)(w p p p + w p p p )d p + p 2 w p w p d p ))fl 



(48) 

In cq.(48) we have 

D p = o^"d v (49) 

where from positive definiteness of — Aw it follows that the momentum depen- 
dent coefficients a M!y satisfy the positivity condition o M o„a'"' > 0. We find the 
equilibrium solution from the requirement 

(£>" + \p 2 P^w v )n E = (50) 

There is a solution (which is the Juttner equilibrium [38]), if 

\ = (3{b-a) (51) 

Then, 

ft E = cxp(-f3wp - jp 2 ), (52) 

where 7 is an arbitrary constant. 

At zero temperature (realized by setting 6 = and A = in eq.(46)) the 
diffusion with a continuous mass distribution (8) can be decomposed into the 
diffusions on the mass shell using the formula 

- J da(p)fd,p 2 P^d,g = J dpdm 2 p{m 2 )5{p 2 - m 2 c 2 )fA%g (53) 
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where 

= (<^m 2 c 2 +^V)^^ + ZpiJLr (54) 

j, fc = 1,2,3. and is the generator of the relativistic diffusion defined by 
Schay[26] and Dudley [27] (see [39] for recent mathematical results). The dif- 
fusion corresponding to b = and A ^ has been derived in [40] (and in a 
covariant form in [41]) from the requirement of the detailed balance leading to 
the Juttncr equilibrium distribution. 

On the basis of the diffusion theory we can develop a statistical description 
of a gas of diffusing particles. We define the current 



= J da(p) P "n (55) 

The current can be expressed by the mean momentum v and the velocity u ( 
u 2 = 1) 

AT" = p(Q)v^ = p{Q)nu^ (56) 
where p(fl) = J da(p)tt, w M = nit 11 and n 2 = v^v^. Then, using 

p^D" + \p 2 P^w u ) = (57) 

we show that 

d*N» = (58) 

It follows that 

N = j dxN (59) 

is a constant (in such a case the probability distribution Q can be normalized). 
Then, let us define the energy-momentum tensor 

9^ = J do(p)p^p v ?L (60) 

From eq.(48) we obtain 

= -{b - 3a)N u - \Q»w v + XQ^w^ - 2bw u w^N^ (61) 

The statistical description can proceed with some thermodynamic notions 
(the definitions follow the ones for mass-shell diffusion of ref.[42], see also [43]). 
We define the relative entropy current 

S k = N ~* I dcr(p)pnn IniN^flNE^ 1 ) (62) 
It can be shown that J dxS 1 ^ > and 

9 i* s k = ~ N_1 J da(p)Qa^d„ In Rd u In R < (63) 
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where R = flfl^ 1 and (denned in eq.(49)) is positive definite . It follows 
that / dxSft is a non-negative function decreasing monotonically to zero (it can 
be interpreted as the entropy of the system plus its heat bath) . The entropy of 
the particle system is 

S = N- 1 J dxj dcr(p)p n\nn (64) 
We define the free energy 

T = $~ X N J d*S° K - N/3- 1 \n(NN E r ) (65) 

and the energy 

W = J dxToo (66) 

It follows from eqs.(64)-(66) that the thermodynamic relation 

(i- 1 S = W-F (67) 

comes out as an identity. The time evolution of each term in eq.(67) is deter- 
mined by eqs.(61) and (63). 

The energy-momentum tensor can be expressed in the form 

9^ = eu»u v + r"" (68) 

where e = p(Q)n 2 and 

= J da{p){p^ - v^){p v - v v )n (69) 

Let 

Ka = Tjua - U v U a (70) 

The divergence of the energy-momentum (68) can be compared with the one 
calculated in eq.(61). As a result we obtain hydrodynamics equations for the 
velocity (relativistic Navier-Stokes, see [34]) 

eurd%u a + h va d*T» v = Xh va (e^w^ - 0/>") - 2bh av w l/ w^ (71) 

If A = b = and t^ u = — 7TE(f2)ft, MI/ (ideal fluids) then in the non-relativistic 
limit eq.(71) coincides with the Euler equation. The terms on the rhs of eq.(71) 
describe a friction which brings the gas of particles to the equilibrium described 
by the state (52)( this is the state of an ideal fluid with the velocity — w M ). 
In order to describe the fluid motion close to the equilibrium we can apply an 
expansion of the solution of eq.(48) in p^d*. In such a case we could determine 
the viscosity by means of the parameters entering the diffusion equation (48) 
which can be calculated in quantum gauge theory at finite temperature. 
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5 Discussion and summary 



There is a basic difficulty in a formulation of a transport theory of quark-gluon 
plasma which starting from quantum gauge theory at finite temperature could 
reach a form useful for calculations. First of all the non-perturbative QCD is 
formulated in the Euclidean framework on the lattice. In such a formulation it is 
impossible to treat quarks dynamically in the real time. Some approximations 
are unavoidable. Our approximation relies on the Wong equation which can 
be derived from formal continuum quantum gauge theory [11] [12]. In this paper 
we have shown that the diffusion results from a model of a particle moving 
in a gauge field in the same way as a non-relativistic diffusion results from an 
evolution in a random electric and electromagnetic fields [17] [18]. The Wong 
equation has the gauge coupling constant as the only parameter. Then, we de- 
rive the diffusion equation which depends on three parameters p,tte and A which 
could be calculated from the correlation functions in the Euclidean framework 
on the lattice (the diffusion constant has been expressed by correlation functions 
of Wilson loops in finite temperature QCD in an earlier paper [16]) . When the 
diffusion equation is applied then it could be directly compared with experimen- 
tal data (as some diffusion equations are [5]) or we could use it for some other 
approximations. We have in particular in mind the hydrodynamic equations 
which could be derived from the diffusion equation. The parameters entering 
the fluid equations could be expressed by the ones of the diffusion equation. 
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